Edinburgh/2002/01 



Semi-leptonic decays of heavy mesons and the 
Isgur-Wise function in quenched lattice QCD 



UKQCD Collaboration 

K.C. Bowler a , G. Douglas^ 1 , R.D. Kenway a , G.N. Lacagnina a , 

CM. Maynard a 

^Department of Physics & Astronomy, University of Edinburgh, Edinburgh EH9 

3JZ, Scotland, UK 



Abstract 

The form factors for the semi-leptonic B —* D and B — * D* decays are evaluated in 
quenched lattice QCD at two different values of the coupling, (3 = 6.0 and 6.2. The 
action and the operators are fully 0(a) non-perturbatively improved. The slope of 
the Isgur-Wise function is evaluated, and found to be p 2 = 0.83l}^l^ 4 (quoted errors 
are statistical and systematic respectively). Ratios of form factors are evaluated and 
compared to experimental determinations. 
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1 Introduction 



The B — > D and B — ► D* semi-leptonic decays are of considerable phenomeno- 
logical interest. They provide ample opportunity for interaction between ex- 
periment and theory. The ideas of Heavy Quark Symmetry (HQS) [1-5] and 
Heavy Quark Effective Theory (HQET) [6] were first developed and applied 
to these decays. A combination of theoretical and experimental data for semi- 
leptonic decays can be used to determine the CKM matrix element \V c b\. HQS 
allows theoretical control of the non-perturbative (NP) aspects of the calcu- 
lation around the infinite quark mass limit. In this study we determine the 
non-perturbative matrix elements directly on the lattice. 
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HQS can be used to constrain the form of the matrix elements that describe 
semi-leptonic decays of heavy-light mesons. In particular, the relevant matrix 
elements are expressed in terms of a set of form factors that contain the non- 
perturbative physics of the decay: 



(D(v')\^\B M )_ = {v + vYh+H + {v _ vYh _ {u>) (1) 



y?m B m D 
(D*(v', e)\c^b\B{v)) 
y/m B m D * 
(D*(v',e)\c^ 5 b\B(v)) 



= te^elv' p vM^) (2) 

=(uJ + iy»h Al (u) + 

-[h A2 (u)v^ + h A3 (u)v'%e*-v) (3) 



where v, v' are the velocities of the initial and the final meson respectively, 
and u — v • v'; e is the polarisation vector of the D* meson. The u variable is 
kinematically constrained to the interval 

l<cu< — S^L . 4 

2mBm D (,) 



In the limit of infinitely heavy bottom and charm quarks, the six form factors 
are related to a universal function known as the Isgur-Wise function [4,5]. 
Away from the heavy quark limit this relation is modified by two kinds of cor- 
rections: perturbative QCD corrections and heavy quark symmetry breaking 
corrections. For large enough heavy quark masses, the relationships between 
the form factors and £(u>) are 

hj(uj) = [aj + (3j(m b ,m c ]u) + ^j(m b ,m c ]u) + 0(l/mlj]^(uj) . (5) 



The Oij terms are constants that fix the behaviour of the form factors in the 
heavy quark limit (a + = ay = a Al = «a 3 = 1; a_ = a A , 2 = 0). The (3j and 7j 
functions account respectively for radiative corrections and power corrections 
proportional to the inverse of the heavy quark mass. The radiative corrections 
are calculable in perturbation theory [7], while the power corrections are non- 
perturbative in nature [8]. At zero recoil (v = v',cu = 1) however, Luke's 
theorem [9] guarantees that 7+ = and j Al = 0, which means that power 
corrections to h+ and h Al are of order 0(l/mjj c ). 

At small recoil, the Isgur-Wise function is modelled by 

eH = i-/)V-i) + %-i) 2 ) (6) 
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where p 2 is called the slope parameter and £(1) = 1, because of current conser- 
vation. Alternative parametrisations of £(o;) are possible, which start to differ 
from (6) at 0((u - if) [10,11]. 



In this paper we present a study of the form factors of the B — > D and 
B — > D* semi-leptonic decays, and the extraction of the Isgur-Wise function 
in the quenched approximation to lattice QCD. The calculations are performed 
at two values of the coupling, [3 = 6.0 and (3 = 6.2, using a non-perturbatively 
improved relativistic Sheikholeslami-Wohlert (SW) [12] fermion action and 
current operators, so that the leading discretisation errors appear at 0(a 2 ) 
rather than 0(a), where a is the lattice spacing. Whilst the improvement 
gives better control over the continuum extrapolation, it does not necessarily 
reduce the size of lattice artifacts at fixed coupling. With only two values of 
the coupling, no continuum extrapolation is attempted. 

The remainder of the paper is organised as follows. In section 2 we give de- 
tails of the calculation including extraction of the form factors. In section 3 
we describe our determination of the vector and axial current renormalisation 
constants. In section 4 we report the results of the extraction of the Isgur- 
Wise function and its slope. Section 5 dwells on the problem of quark mass 
dependence, while section 6 describes a calculation of ratios of form factors 
and its comparison with experimental determinations. An analysis of the sys- 
tematic sources of uncertainty and a comparison with other results is detailed 
in section 7. 



2 Details of the calculation 

2.1 Improvement of action, currents and masses 
The improved action used in this work has the form 



where Sw is the Wilson action. The improvement programme also requires 
the improvement of current operators; in particular, the improved vector and 
axial currents are 



(7) 



X 



V*(x) = V^x) + ac v d u T^(x) 
A I li (x)=A li (x) + ae A d li P(x) 



(8) 
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where 



Ap(x)='$(x)'y li 'y 5 i/>(x) 
P(x) =4>(x)>y 5 ^(x) 
T^ v (x) =4>(x)icr llu il)(x) 



and is the symmetric lattice derivative in the fi direction, defined in terms 
of the lattice spacing a by 

dj(x) = ^[f(x + ajl)-f(x-afi)} . (9) 



The generic renormalisation of improved current operators is performed as 
follows (J = A,V): 

J R = Zj(g 2 )(l + bj(g 2 )am q )J l (10) 



where Zj is calculated in a mass-independent renormalisation scheme. It is 
conventional to define effective renormalisation constants: 

Zf = Zj{\ + bjam q ) . (11) 



The bare quark mass, in terms of hopping parameters, is equal to 

2 \ K q ^crit / 



where K crit is the value of the hopping parameter at which the bare quark mass 
vanishes. The improved, renormalisation group invariant (RGI) quark mass is 
defined in the following way 

am q = Z m am q {l + b m am q ) . (13) 



A discussion of the different determinations of the improvement coefficients 
can be found in [13]. We use the NP value of cgw determined by the ALPHA 
collaboration [14,15] and the improvement coefficients for the current opera- 
tors are taken from Bhattacharya et al. [16-18], except for Z rn and b m . The 
additive b m coefficient is evaluated at 1-loop in perturbation theory [19] with 
the coupling "boosted" by the mean link, g 2 — > g 2 /uQ, and Z m is determined 
non-perturbatively by the ALPHA collaboration [20]. 
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Table 1 

Simulation details. 



= 6.2 = 6.0 

Volume 24 4 x 48 16 3 x 48 

iV C onfig S 216 305 

c SW 1-614 1.769 

a- 1 (r )(MeV) 2.913 2.123 

K hcavy 0.120,0.1233,0.1266,0.1299 0.1123,0.1173,0.1223,0.1273 

/Slight 0.1346,0.1351 0.13344,0.13417 



2.2 Simulation details 



The lattice spacing is set using the Sommer scale r [21,22]. The matrix el- 
ements are extracted from combinations of quark propagators corresponding 
to four values of the heavy quark masses (m^eavy — m c) and two values of 
the mass of the light (passive) quark (mii g ht — iris). Simulation details are 
summarised in Table 1. The gauge configurations that are used in this calcu- 
lation were generated with a combination of the over-relaxed [23,24] and the 
Cabibbo-Marinari [25] algorithms with periodic boundary conditions. The 
light quark propagators are smeared with the fuzzing technique [26], while 
heavy quark propagators are smeared using a gauge invariant technique [27]. 
The statistical errors are estimated using the bootstrap method [28]. 



2. 3 Extraction of the form factors 



The matrix elements relevant to the semi-leptonic decays are extracted from 
fits to two and three-point correlation functions. The general form of the latter 
is as follows 

c^j(p A ,t x ;p B ,t y ) = x;e- i ^-^(o|r[nfl(x)j"t(y)ni(o)]|o> (14) 



where q '= p B — Pa and Q x is the operator that creates the state X. J M is the 
weak current. We have computed the correlation function using the "standard 
source" method [29,30]. The time-slice t x is fixed to t x = 28, i.e. slightly off the 
centre of the lattice, in order to allow for an analysis of the different systematic 
effects of the two sides of the lattice. The extended heavy quark propagators 
were computed at four values of the hopping parameter. For the active quark 
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propagators, which are part of the current, only two values of the hopping 
parameter were used at (5 — 6.2 due to disk space constraints and all four 
were used at (5 — 6.0. 

In the Euclidean time formulation, when the time separation of the operators 
is large enough, i.e. when t y and t x — t y are large, one has 



where Zj = (0|Pj|Pj(pj)). The Zi constants and the energies Pj = \J\pi\ 2 + rnf 
are extracted from fits of two-point functions (detailed descriptions of the 
methods employed can be found in [13,31]). The {Z iy Pj} parameters are then 
used to cancel the asymptotic time dependence from the three-point function 
and to extract the desired matrix element from a fit to a plateau. For each 
value of uj where 

u = m% + ml w -<? 
2m B m D( .) 



an average is performed of all the kinematic and Lorentz channels for which 
the matrix elements have the same value. An example matrix element is shown 
in figure 1. All distinct matrix elements are fitted simultaneously to extract 
the form factors. We examine five kinematic channels, each of them being 
specified by the values of \pa\, \q\- These are summarised in Table 2. 

Since the heavy quarks used in this study have masses around that of the 
charm quark, we are in effect extracting the Isgur-Wise function from aD^ 
P/*) decay. Provided that there is little residual mass dependence beyond 
leading order power and radiative corrections in the form factors, then the 
machinery of HQS allows us to compare to experimental B — > I)W decays. 

Table 2 

Spatial moduli of momenta in the different kinematic channels, in lattice units of 
2ir/(aL). 
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Wa\ 













1 





i 


2 


1 


i 


3 


1 





4 


1 


V2 



6 



0.003 



0.002 



0.001 



1 I 1 I 1 I 

i 



i 1 i 1 i 



iiiiiiii 
****** 



_i i i i i_ 



i.i.i 



i.i.i 



i.i.i 



2 4 6 



10 12 14 16 18 20 22 24 26 28 



0.02 



S-H 

cL 0.015 

s 

CD 



0.01 



f 1 I 1 I 1 I 



i i i i i i i i i i i i i i i i r 



I 1 I 1 I 



i . i . i . i . i . i . i . i 



i . i 



j_ 



2 4 6 8 10 12 14 16 18 20 22 24 26 28 

time 

Fig. 1. The matrix element (P b (p b )\V^\P a (pa)) with (3 = 6.2, k h = k h > = 0.1200, 
kl = 0.1346, \pa\ = and \pb\ = 1- The data shows the ratio of the three-point 
function over the fitted two-point parameters. The circles are < t y < 27 (fore side) 
and the diamonds are < (T — t y ) < 19 (back side). 

3 Current renormalisation 



3.1 The vector current 



Observing that the degenerate mass vector current is conserved in the con- 
tinuum, we consider the effective matching coefficient , obtained from the 
forward matrix element of the current between degenerate heavy-light mesons 
at rest. We then have 



(P(0) \V \P{0) ) cont = 2M = Zf(P{6) |Vtf|P(0) ) 



latt 



(17) 



In particular, taking the ratio of a heavy-light two-point function, at the ex- 
tension time-slice, with the forward degenerate matrix element gives, 



C 2p t(p = Q;t = t a 



2„-Mt x 



Z 2 e 



C^ B (PA=ph = 0;t x ,t y ) 2M 
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Table 3 

Zy 1 and quark mass. arriQ is the bare quark mass appearing in the current and am q 
is the bare quark mass of the light, spectator quark. 





1 


(Win 








am — 0339 


nm — n m Q5 


6.2 


0.485 
0.268 


1 ^1 99+ 3 

1.087ll| 


1 ^1 1 1 + 3 
l.Olll ^ 

1.08731^ 






om, = 0.0502 


am, = 0.0298 




0.756 


1.5963l| 


1.59461^ 


6.0 


0.566 


1.3876l| 


1.38431^ 


0.392 


1.20581^ 


1.2044l| 




0.231 


1.0333±| 


1.0333±| 


(2M) 2 



(P(0) 1^1^(0) )latt^ e -Mt HZe -M(t x -t y ) 
2M 

~ (P(0) 1^1^(0) ) latt 

_ VCff 
— Zy . 

The results as a function of bare quark mass, as it appears in equation (11), 
are given in Table 3. Also shown in the table is the spectator quark mass. 
We find that the current renormalisation does not depend significantly on the 
spectator quark mass. The results for /3 = 6.2 have already been reported in 
[32] . We note a small discrepancy that is accounted for by the following: in [32] 
we used the fitted two-point function values of the mass and overlap factor, 
whereas here we use the two-point function itself to determine Z^. This is 
expected to reduce the statistical errors as the ratio is independent of t y . 



3.2 The axial current 

The axial current is not conserved, and so we rely on HQS to determine the 
effective renormalisation constant. The form factor h,A x at zero recoil is equal 
to unity, up to symmetry breaking corrections, which are reduced to 0(1/i71q) 
by Luke's theorem. This form factor can be measured on the lattice, and is 
related to the effective renormalisation by, 

Zftig(l) = l+(3 Al (l)+0(l/m 2 Q ) . (19) 
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(3 = 6.2 (3 = 6.0 




Fig. 2. Zy versus bare quark mass at both values of f3. am c is the charm quark 
mass in lattice units, as extracted from a lattice determination of the heavy- light 
pseudoscalar meson mass. The range indicates the variation in the mass due to the 
different choices of the scale-fixing quantity. 

To measure the form factor, we do not rely on degenerate initial and final 
states only. When they are not degenerate, the renormalisation depends on 
the average mass of the quarks in the current. Our results are shown in Table 
4. Entries that have the repeated average quark mass are taken from matrix 
elements with the initial and final states interchanged. 

3.3 Comparison with other determinations 

We compare our measurement of the effective renormalisation to equation 
(11) using values of Zj and bj determined non-perturbatively by the ALPHA 
collaboration [33] and Bhattacharya et al. [16-18]. This comparison for the 
vector current is shown in Figure 2 for the heaviest spectator quark. The 
agreement is quite striking, especially considering that the ALPHA collabo- 
ration determined the coefficients at near zero quark mass and Bhattacharya 
et al for light quarks (strange scale). Indeed, Bhattacharya et al. have already 
compared their results at (3 = 6.2 to those in [32]. 

For the axial current, there is no NP determination of b A by the ALPHA 
collaboration. To facilitate a comparison, we use the one-loop value of b A [19] 

b A = 1 + 0.1522 g 2 (20) 



9 



Table 4 

and average heavy quark mass. 



g K) zf 

am q = 0.0332 am q = 0.0195 



6.2 



0.485 


1.29+1 


1.2711 


0.429 


1.2611 


1.251* 


0.377 


1.24+| 


1.221* 


0.377 


1.23+J 


1.23+ 4 


0.326 


1.20±| 


1.19l 3 3 


0.321 


1.18l| 


1.171* 


0.268 


i-io±| 


1.0911 


0.218 


l.ioli 


1.0911 



am q = 0.0502 am q = 0.0298 



6.0 



0.756 


1.55+1 


1.521? 


0.661 


1-4511 


1.43l| 


0.661 


1.48+J 


1.46l^ 


0.573 


1.4011 


1.38l| 


0.573 


1.431! 


1.4215 


0.566 


1.371 4 


1.361! 


0.493 


1.33+ 4 


1.311! 


0.493 


1.4111 


1.4115 


0.479 


1.321 4 


1.311! 


0.479 


1.33+ 4 


1.321! 


0.398 


1.271 4 


1.25+ 4 


0.398 


1.201 4 


1.321! 


0.391 


1.22+ 4 


1.211! 


0.311 


1.191 4 


1.18+ 4 


0.311 


1.231 4 


1.231! 


0.231 


1.0311 


l.Oll 4 



with the "boosted" coupling, g 2 = go/u^, where the mean link, uq, is evaluated 
from the 4th root of the average plaquette. This is shown in Figure 3. We 
also show the values of evaluated in boosted perturbation theory (BPT). 
Again, our data compares very well to other determinations. 
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(3 = 6.2 (3 = 6.0 




0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

<am^> <am_> 



Fig. 3. against average bare quark mass. 

We don't use our results to determine bj and Zj\ however, the comparison 
suggests that lattice artifacts are not spoiling the renormalisation. 



4 Extraction of the Isgur-Wise function 



The form factors that are protected from power corrections by Luke's theorem 
are also used to extract the Isgur-Wise function: 



™°TT%b> (21) 



The power corrections for h + and Ha 1 have been estimated and found to be 
consistent with zero in the range 1.0 < uj < 1.2, and therefore are neglected. 
It is not possible to extract the Isgur-Wise function from the form factors for 
which otj = as they differ from zero only by power corrections. In principle 
it would be possible to measure the power corrections in this way. However, 
the poor quality of the data for the sub-leading form factors prevents this. 
Similarly, the hy form factor could not be used because its power correction 
is found to be large, as expected [7]. 
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Fig. 4. The Isgur-Wise function from the B — > D and B — > L>* semi-leptonic decays, 
with Kii g ht = 0.1351, including single-decay and simultaneous fits. 



Table 5 

Results of the linear fits to the Isgur-Wise function, at [3 = 6.2. Quoted errors are 
statistical. 





B - 


-> D 


B - 


> D* 


Sim. 


fit 


flight 


0.1346 


0.1351 


0.1346 


0.1351 


0.1346 


0.1351 


P 2 


0.97±* 4 


0M+\l 


1.431" 


0.93^ 9 7 


1.06±iJ 


0.83±if 



The Isgur-Wise function is fitted to the linear model (6), at both values of the 
coupling. Simultaneous fits of the Isgur-Wise function from the two distinct 
decays are also performed; the assumption is made that the effect of mixing 
different systematic errors (different decays) is no larger than the effect that 
one gets by fitting simultaneously data from correlators with different mo- 
menta. Results are summarised in Tables 5,6 and plotted in Figure 4. The 
determination of the slope of the Isgur-Wise function at (3 — 6.2 from the 
simultaneous fit of the two form factors, at the lightest value of the passive 
quark mass, is taken as the central value. 
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Table 6 

Results of the linear fits to the Isgur-Wise function, at (3 = 6.0. Quoted errors are 
statistical. 





B - 


■+ D 


B 


D* 


Sim. 


fit 


flight 


0.13344 


0.13417 


0.13344 


0.13417 


0.13344 


0.13417 


P 2 


0.94±?{j 


0.9ll| 


0.681^ 


0.88i^ 


0.88111 


0.90l| 



5 Quark mass dependence 

Close to the heavy quark limit (mb,m c — > oo), at each fixed value of uj, 
the Isgur-Wise function does not depend on the mass of the heavy quarks. 
However, the finite size of the simulated masses can introduce a quark mass 
dependence, which has to be examined. 

5.1 Light quark mass dependence 

In the dataset used in this work, only two values of the passive quark mass 
are available at each value of the coupling, making light quark mass extrapo- 
lations impossible. The values of the form factors with the two different light 
masses are compared and found to be statistically consistent with each other. 
However, the global fits of the Isgur-Wise function with the two different pas- 
sive quarks yield slightly different, if still statistically compatible results. This 
effect is taken into account as a systematic uncertainty on the extracted slope 
of the Isgur-Wise function. 

Figures 5 and 6 show a comparison of the fits of the Isgur-Wise function with 
the two different values of the light quark mass, at f3 — 6.2, for both decays. 

5.2 Heavy quark mass dependence 

In this work the radiative corrections are applied, and the Isgur-Wise func- 
tion is evaluated in a range in uj in which the lattice estimates of the power 
corrections for h + and Iia 1 are found to be statistically consistent with zero. 
The Isgur-Wise function is extracted from the form factors by dividing off the 
radiative corrections. This procedure should remove all residual quark mass 
dependence at fixed uj, if we are sufficiently close to the heavy quark limit. To 
study any residual heavy quark mass dependence, the Isgur-Wise function is 
extracted from all the momentum channels, holding all the quark masses fixed. 
For each heavy quark mass combination, the Isgur-Wise function is fitted, us- 
ing all the kinematic channels, and its slope is extracted. Then, for each value 
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Fig. 5. The Isgur-Wise function from the form factor, for both values of the 
light quark mass, at (5 = 6.2. 
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Fig. 6. The Isgur-Wise function from the /i+ form factor, for both values of the light 
quark mass, at (3 = 6.2. 
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Fig. 7. Fits of the Isgur-Wise function at fixed values of the active quark propagator 
mass against the mass of the extended quark propagator, with Kught = 0.1346, at 
(3 = 6.2. Quoted quark masses are improved and renormalised according to (13). 

of the extended quark mass, the result of the fit is plotted against the value 
of the active quark mass. The results, for both values of (3, at fixed values of 
the light quark mass, are shown in Figures 7 and 8. As one can see, there is 
no statistically significant heavy quark mass dependence. It should be noted 
that fixing the quark masses considerably reduces the number of points that 
are available to each of the fits: consequently, some of the fits are quite poor. 



6 Ratios of form factors 

The two following ratios of form factors are also calculated: 



These two ratios would be equal to one in the absence of symmetry breaking 
corrections. Figures 9 and 10 show that this work's determinations are in 
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Fig. 8. Fits of the Isgur-Wise function at fixed values of the active quark propagator 
mass against the mass of the extended quark propagator, with Ku g ht = 0.13344, at 
(3 = 6.0. Quoted quark masses are improved and renormalised according to (13). 

agreement with the experimental determinations by the CLEO collaboration 
[34], that quote 



#1 = 1.18(32) 
#2 = 0.71(23) . 



(25) 



In the case of R\ there is some evidence of a systematic deviation between our 
results and experiment. However, it should be noted that this form factor has 
a strong dependence on the improvement coefficient of the vector current, cy, 
which is poorly known even at (3 = 6.2. 



7 Systematic uncertainties and comparison with other results 



7.1 Systematic uncertainties 



Several sources of systematic uncertainty are considered. Discretisation effects 
are estimated by comparing the results at the two values of (3. The results at 
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Fig. 9. The ratio R%. The range of the experimental determination is shown by the 
dashed lines. 

(3 = 6.2 are taken as our central values. At (5 — 6.0, the extraction of the form 
factors is affected by considerably larger statistical errors. This is reflected in 
the large statistical errors on the slope of the Isgur-Wise function. 

The systematic uncertainty induced by the scale fixing quantity is evaluated by 
comparing fits of the Isgur-Wise function using r and m p to set the scale. In 
all cases, it is found to be an effect smaller than 1%. This is mainly due to the 
fact that the radiative corrections depend only logarithmically on the lattice 
spacing (through quark masses). Furthermore, the effect of the definition of 
the RGI quark mass and the values of the coefficients Z m and b m only enters 
through the radiative corrections. 

The quark masses used in this calculation are around the charm quark mass. 
It is necessary to check for any residual heavy quark mass dependence beyond 
radiative and leading order power corrections. This has been examined and 
no statistically significant trend has been observed (see section 5.2). No nu- 
merical estimate for this is included in our systematic error. The light quark 
masses are around the strange quark mass, and thus the Isgur-Wise function 
extracted is £ s rather than £. We examine the light quark mass dependence 
by comparing £ measured at the two values of quark mass. An extrapolation 
is not attempted, but the lightest mass is taken as a central value and the 
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Fig. 10. The ratio i?2- The range of the experimental determination is shown by the 
dashed lines. The data points are scattered in a wide range because they correspond 
to different kinematic channels, with different systematic effects. 

difference systematic uncertainty. 

The statistical errors on are much larger than those on h + and so the 
latter form factor dominates the simultaneous fit. Fitting to either form factor 
on its own produces the same result within statistical error. In all cases, the 
data satisfy the HQS condition £(1) = 1 to such an extent that imposing this 
condition on the fit has no sizeable effect. 

It is rather difficult to quantify the effects of quenching other than by perform- 
ing the calculation in full QCD. There is some evidence that quenching could 
have a sizeable effect in heavy-light decay constants [35], around 10 — 15%. 
However, the latter are dimensionful quantities and the scale setting ambiguity 
in the quenched approximation is sizeable. It is possible that large unquenching 
effects in the decay constant are at least partly due to the scale setting ambi- 
guity [36]. The Isgur-Wise function is a dimensionless quantity and we have 
seen that it is almost independent of the choice of the scale fixing quantity. 
There is some evidence that quenching effects on the ratio are mild [37]. 
It is reasonable to assume that quenching effects will not be larger than any 
other systematic effect in this calculation. 

Our systematic uncertainties are then the difference of £ obtained with the 
two values of the light quark mass (+23%), with the two values of the cou- 
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pling (+8%) and with the quantity used to set the scale (—1%). Positive and 
negative systematic uncertainties have been separately added in quadrature. 



Therefore, our final result is 



p 2 = 0.83^t 24 . 



(26) 



7.2 Comparison with other results 

Other recent lattice results for heavy-light to heavy-light decays have concen- 
trated on computing the value of the physical form factor at zero recoil to 
extract \V c b\. They use ratios of matrix elements to cancel uncertainties [38] 
and further calculate the HQS breaking corrections [39]. Preliminary UKQCD 
results on this dataset for h + [40] and h J y 1 [41] are consistent with each other 
and this work. Previous UKQCD determinations [42-44] give p 2 = l^lj^i 
with the spectator quark mass around strange, and p 2 = O.Oi 2 ^ when ex- 
trapolated to the chiral limit. Hashimoto and Matsufuru [45] compute £ using 
lattice HQET, but with an incomplete set of the 1/itiq corrections, and quote 



We also compare our result for the slope of the Isgur-Wise function to ex- 
periment. Both CLEO and Belle have recent measurements of B — > D* de- 
cays [46,47] and quote 



Our result is significantly lower than the experimental ones, but consistent 
with other quenched lattice determinations. However, the experiments mea- 
sure the slope of J-{uS), the physical form factor, which is equal to the Isgur- 
Wise function only in the infinite quark mass limit. Applying radiative cor- 
rections t]a to our results would increase the slope by around 5%. Considering 
the size of our uncertainties, this effect is small and we neglect it for this 
comparison. 
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